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X-ray Variability Coherence in the Time Domain 

Xiao Che\ Cong-Zhan Liu^ and Ti-Pei Li^-^-^ 
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I A technique for calculation of variability coherence at different timescales 

pL( i performed directly in the time domain is introduced. Simulations are made to 

Q\ I compare the coherence spectrum derived by the time domain technique with the 

^ I coherence function by Fourier analysis. The results indicate that in comparison 

>• I with the Fourier analysis the time domain technique is more powerful for revealing 

^ I signal coherence in noisy data. We apply the time domain technique to the real 

' data of the black hole binaries Cygnus X-1 and GX 339-4 and compare the results 

Q ' with their Fourier coherence spectra. 

i> : 
o ■ 
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I 1. Introduction 

. 

^ ^ ■ The complex aperiodic variability of high-energy emission shown in different time scales 

^ ' is a common character for X-ray binaries and super massive black holes. Temporal anal- 
■ ysis is an important approach to study undergoing physics in objects from observed light 
curves. The spectral lag between variation signals in different energy bands provides useful 
information on their producing and propagating processes in celestial objects. Besides the 
time lag, the spectral coherence, a measure of the degree of linear correlation between two 
light curves in different energy bands, is another useful quantity to provide constraints on 
models of physical processes in observed objects. For studying variability coherence in the 
frequency domain, the coherency function - the normalized average amplitude of the cross 
power density as a function of Fourier frequency - can be used. Vaughan and Nowak (1997) 
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have presented a procedure for estimating the coherence function in the presence of counting 
noise, which has been apphed to studying black hole candidates (e.g. Nowak & Vaughan 
1996; Vaughan & Nowak 1997; Cui, Zhang & Chen 2000; Nowak et al. 2001; Nowak, Wilms 
& Dove 2002; Ji et al. 2003). 

The time domain quantity corresponding to the coherence in the frequency domain is 
the coefficient of linear correlation between two light curves in different energy bands. The 
correlation coefficient has not been extensively applied in temporal analysis as the coherence 
function in the frequency domain. One of the reasons is that the variability is caused by 
various physical processes at different timescales, a single value of correlation coefficient from 
two observed light curves is not enough to provide measures of association between complex 
variables. What needed is a distribution of correlation coefficient over timescales. Here we 
present a procedure for deriving the coherence spectrum in the time domain, or a timescale 
spectrum of the correlation coefficient for a complex process in the presence of background 
noise and compare the coherence spectrum in the time domain with the Fourier frequency- 
dependent function of coherence by simulation in §2. In our simulation, the result from the 
time domain technique is consistent with the coherence function obtained from the Fourier 
analysis when the signal to noise ratio is high. In general, the time domain technique can 
provide a more accurate and stable result than the Fourier analysis. We apply the time 
domain technique to the black hole X-ray binaries Cygnus X-1 and GX 339-1 and present 
the results in §3 and give discussion on timelag effect in §4. 



A proper statistic for measuring the degree of linear correlation between variations in 
two A^-point signal counting series Xs{i) and ys{i) is their correlation coefficient 



For two observed lightcurves x{i) — Xs{i) + Xn{i) and y{i) = Usii) + Unii) with independent 
noise x„ and y„ respectively, we have 



2. Timescale Spectrum of Correlation Coefficient 



a{xs)cr{ys) 




^{x{i) - x){y{i) -y)^ ^{xs{i) - Xs){ys{i) - ys), 



and 



a'^(x) = a'^ixs) + (T^{xn), a'^iy) = cr'^iys) + cr^iyn), 
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then the signal's correlation coefficient can be calculated as 

^ ^ Y^tii^sji) -^s)iysii) -ys)/N ^ T.^=Ax{i)-x){y{i)-y)/N 

<^s)c^{ys) ^a\x) - a\xn)sja^{y) - a^yn) ' 

It should be pointed out that one can not take two so called "background-subtracted lightcurves" 
as signal's lightcurves Xg and yg to calculate the correlation coefficient simply by Eq. (1), be- 
cause background lightcurves given by a database are usually smoothed curves not including 
the fluctuation of real noise at shorter timescales. 

We apply a timescale analysis technique (Li, Feng & Chen 1999; Li 2001; Li et al. 2004) 
to calculate correlation coefficients for different timescales from two observed lightcurves. 
From an original lightcurve x{j; St) {j = 1, 2, ■ ■ ■) with a time resolution 6t, we can construct 
M different lightcurves with the same larger time step At = M6t by combining M successive 
bins 

iM+m-l 

Xm{i]At)= Yl <r,5t) , (m = l,---,M), (3) 

j=(j-l)M+m 

where the combination starts from the mth bin of the original lightcurve. From two simul- 
taneously observed counting series x{i] 5t) and y{i] 6t), the correlation coefficient r{At) for 
a certain timescale At — M^tSt can be calculated by 

The total observation period is divided into L segments. The correlation coefficient rk{At) 
is calculated by Eq. (4) for a pair of lightcurves on each time segment k separately. The 
average f(At) = J2k=i ^k{At)l L and its standard deviation aij) are defined as the coherence 
coefficient for the timescale At and its error. 

To check the procedure described above and compare it with the Fourier coherence 
function we produce two 1000 s signal time series, Xg and ys, with time bin 5t — 1 vas 
sampled from a two-dimensional normal distribution 

p{xs,ys) — exp 

27ra{xs)a{ys)^/l - 

where (xg), cr{xg), {yg),cr{yg) are means and standard deviations of Xg and yg respectively, and 
p is the correlation coefficient between Xg and yg. By giving these parameters in Eq. (5), we 
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are able to produce two signal time series with a certain correlation coefficient. To simplify, 
we specify (xg) = (Vs) = 49, a{xs) = (j{ys) = 7. The correlation coefficient p in Eq. (5) is the 
theoretical value of coherence, and it is used to compare with results of time domain coherence 
and Fourier coherence function. We add independent Poisson noise a;„ and ?/„ on the signal 
series Xg and Us respectively to get simulated lightcurves x = Xs + Xn and y = i/s + Vn- 
Different signal to noise ratio (SNR) could be set by adjusting the mean of independent 
Poisson noise. The simulated hghtcurve is divided into ten segments, each has a length of 
100 s. For each segment, coherences at different timescales are calculated by Eq. (4) and 
Fourier coherence coefficients at different frequencies are also computed with the procedure 
described by Vaughan and Nowak (1997). Finally, the average and its standard deviation are 
computed from the ten time domain coherences and Fourier coherence coefficients separately. 

Figure 1 shows coherence spectra from simulated light curves with p — 0.8 and SNR 
= 5, 2 and 1 separately, where Fourier coherence function in the frequency domain are 
transformed into the time domain by timescale = 1/ frequency . From Fig. 1 we can see 
that, in comparing with Fourier coherence functions, time domain coherence distributions 
are more consistent and accurate. For the case of low SNR, Fourier coherence functions 
fluctuate seriously and have a large error bar, while the time domain procedure has a better 
capability to reveal a given coherence, especially at short timescales. Figure 2 show the 
results for a certain SNR = 1 and different p — 0.9, 0.8 and 0.6 respectively, indicating that 
the time domain procedure is more powerful to reveal lower measured coherence. 

The test and comparison arc extended to different signal distributions other than the bi- 
normal one. The joint probability density function of binormal distribution can be factorized 
into two factors 

p{x,y) ^p{x)p{y\x) , (6) 

where 



p{x) — 



27r(7i 



and 



p{y\x) = 



■ exp 



exp 



1 f X — {x) 

'2 



1 (y-{y)- p^c^: 



(7) 



(8) 



The first factor p{x) of Eq. (6) is a onc-dimcnsional Gaussian density function of the variable 
X which is independent of y. The second factor p{ii\x) is also a Gaussian density but with 
average y + p^^cFy, standard error Oy^Jl — and the correlation coefficient p between 
x and y, which can be seen as a conditional density of y given x. With Eq. (6) we can 
produce a pair of time series with a given correlation coefficient p with a density function 
p{x) different to Gaussian by sampling x from the density p{x) and y from the density of 
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Eq. (8). The exponential distribution, uniform distribution and Poisson distribution are 

adopted as p{x) to produce correlated pairs of time series separately, and then compute 
their coherence distributions in the time domain and their Fourier correlation functions. All 
results arc similar with Fig. 1 and Fig. 2, showing that the time domain procedure is an 
effective approach to coherence analysis. 

3. X-ray Coherence of Black Hole Binciries 

Here wc apply the time domain technique of coherence analysis to two black hole candi- 
dates, Cyg X-1 in the low state and GX 339-4 in the low and high states. We use observational 
data from PC A/ RXTE. PCA/ RXTE observation has dead time 10/iS much smaller than 
time resolution 1ms of original lightcurvcs wc use in this paper. And for Cyg X-1, dead-time 
affection is about 1% (Maccarone and Coppi 2000), so it could be ignored. 

Variabilities in lightcurves observed for an astronomical object are usually caused by 
various processes at different timcscales with different coherences. The coherence coefficient 
r{ At) calculated by Eq. (4) reflects not only the variation property on the timescale At, but 
is also affected by that on larger ones up to the total time period used in the calculation. 
To limit the effect of larger timescales we divide the total observation period into segments, 
each has a duration just lOAt. The resultant correlation coefficient should be then related 
to the limited time scale region of At — 10 At. 

In our calculation for Cyg X-1 and GX 339-4 data, for a given time bin At the studied 
observation period is divided into segments with a duration of 10 At each and the coherence 
coefficient is calculated for each segment by using Eq. (4) . The all segments are divided into 
ten groups with the segment numbers in each group being nearly equal to each other. We 
calculate 10 avarage coherence coefficients rk{At) for each group k {k = 1, . . . , 10) respec- 
tively and take their average r{At) = X]fc=i standard deviation a[r{At)] 
as the final result of the observation. The Fourier coherence coefficients and their average 
and the standard deviation of the average arc also calculated from the 10 time groups. The 
results are shown in Figure 3. For both Cyg X-1 and GX 339-4 in the low state shown in the 
left and middle panels of Fig. 3, both Fourier coherence and time domain coherence show 
similar trend of change from near-perfect coherence at large time scale to low coherence at 
small time scale. Two other black hole candidates, XTE J1550-564 (Cui, Zhang and Chen 
2000) and GRS 1915-1-105 (Ji et al. 2003), have Fourier coherence function with near- unity 
coherence at large time scale and low coherence at small time scale. We also apply the time 
domain technique to these sources and get similar results. As the coherence is much more 
easily to be lost than maintained, only when one linear physical process is involved could 
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coherence be preserved (Vaughan & Nowak 1997). Any nonlinear process or more than one 
independent linear processes involved could degrade coherence. The near-perfect coherence 
at large time scale may provide a strong evidence that there is a powerful physical process 
in these sources with a large characteristic timescale dominates the emission of photons in 
the studied energy bands. 

It is a remarkable feature that all the values of time domain coherence are very low in 
the whole studied timescale range for GX 339-4 in the high state, shown in the right panel 
of Fig. 3, indicating that there is only little relationship between variabilities of 2-5 keV and 
5-7 keV emissions of GX339-4 in its high state. 

To see whether the relationship between variabilities of two light curves coincides with 
the time domain coherence computed from them, we simultaneously plot two observed light 
curves of Cyg X-1 in two energy bands for different time scales at the same time period, 
shown in Figure 4. From Fig. 4 one can see that the time domain coherence calculated by 
Eq. (4) can reflect the real strength of the hnear relationship between two hght curves. 



4. Discussion 

The Fourier coherence function r(/) of two lightcurves x and y is defined in the space 
of Fourier frequency 

.(/) = ^ l<^'/')l . (9) 

vwumwuwi 

Where, X{f) and Y{f) are the Fourier transforms of x and y respectively, the cross spectrum 
C(/) — X*{f)Y{f), and angle brackets denote an average over the used segments of hght 
curve. If the ratio of two Fourier transforms at a frequency /, H{f) — X{f)/Y{f), is the 
same for all segments of the two processes, or, equivalently 

y{t) = Jh{t- T)x{T)dr , (10) 

then r(/) = 1, the processes are said to be coherent at frequency /. But the unit coherence 
in the time domain, r{At) = 1, indicates a linear correlation between the variabilities of 
intensity in the two energy bands on the timescale At 

y{i;At)^hx{i;At) , {t = 1, . . . , Nt) , (11) 

where his a constant during the observation period T — NrAt, which is a stronger constraint 
than Eq. (10) from an unit coherence, r(/) = 1, in the frequency domain. It is an advantage 
for Fourier coherence function that it can obtain coherence at all frequencies simultaneously 
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and exactly, because amplitude at all frequencies could be obtained independently from 
Fourier transformation of original lightcurve. Improved Fourier coherence function (Vaughan 
& Nowak 1997) can calculate coherence of noisy lightcurves and uncertainty as well. While 
the coherence coefficient calculated in the time domain by Eq. (4) can be confined in only 
a timescale region, At — NAt with being the number of time bins in a calculated time 
segment, and the error of the average coherence coefficient is estimated by the statistical 
deviation in coherence coefficients of different time segments. 

Spectral time lags arc usually existed in observed processes. To calculate the time 
domain coherence r{At), one should firstly derive their timescale spectrum of time lag, then 
make time lag correction for light curves at the timescale At. But, in practice, the coherence 
derived by the time domain technique is not sensitive to the spectral time lag in analyzing 
real observational data. From lightcurves of Cyg X-1 in its low state in energy band 2-5 
keV and 5-13 keV observed by PCA/RXTE, we use the modified cross-correlate function 
(MCCF) (Li, Feng & Chen 1999; Li ct al. 2004) to calculate the spectral time lag, the result 
is similar to what obtained from Fourier analysis (Nowak et al 1999). Then we calculate 
time domain coherences between the two lightcurves with and without time lag correlation 
separately. The result shown in Fig. 5 does not indicate any essential effects of the spectral 
lag on the time domain coherence. To see how seriously time lag could affect results of time 
domain coherence, we make a modification to the original observed lightcurves to produce an 
artificial deviation: the band 5-13 keV lags behind the band 2-5 keV by 10 times of time lag 
measured from the observed fightcurves by MCCF, then calculate time domain coherences 
between the modified fightcurves. As we can see in Fig. 5, only a litter deviation can be 
found between the result from the original lightcurves and the modified ones. We also test 
other data used in this paper to see whether Fig. 5 is only an exception, the results turns 
out to be that although the spectral time lag exists in most data, it has just little infiuence 
to the coherence obtained by the time domain technique. 

To understand reasons of insensitivity of time domain coherence to spectral timelag 
in observed X-ray data, we make simulation studies. Two 1000 s lightcurves with 1 ms 
time bin and a given correlation coefficient p = 0.8 (or expected coherence = 0.64) are 
produced by sampling from the two-dimensional normal distribution Eq. (5) described in 
§2. One of the lightcurves is then artificially shifted by 10 ms. We use MCCF to derive 
timescale distribution of time lags and calculate time domain coherences between the two 
lightcurves with and without timelag correction. The results are shown in the left panel of 
Figure 6. We can see that, though the effect of spectral timelag can not be neglected in 
general, the difference of coherences obtained with and without timelag correction is small 
at large timescale. For deriving time domain coherences at large timescales, we need to 
reconstruct original lightcurves with a larger time bin (see the details in §2), which impairs 
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the effect of a timelag which is much smaller than the time bin. Spectral timelags measured 
at a timescale are usually much shorter than the timescale itself. For example, the X-ray 
time lag of Cyg X-1 between 2-5 kcV and 5-13 kcV at the timescale of 1 s (or frequency of 1 
Hz) is only about 10 ms. From the left panel of Fig. 6 we can see that the effect of timelag 
on time domain coherence is small when timelag is much smaller than time bin. 

The correlation existed in successive time bins in observed lightcurves may be another 
reason of the insensitivity of time domain coherence to spectral timelag in observed X-ray 
data. The right panel of Fig. 6 shows results from two simulated lightcurves based on shot 
model with 10 ms timelag. We produce a 1000 s lightcurve with 1 ms time bin consisted of 
random shots similar to what observed from Cyg X-1 (Feng, Li & Chen 1999). The front of 
the shot is represented by a single exponential and the time constant of both ascending and 
descending front is 70 ms. The interval between two successive shots follows an exponential 
distribution with an average of 450 ms. The above produced hghtcurve is shifted by 10 
ms to create the second one. The right panel of Fig. 6 shows coherence spectra derived 
by the two lightcurves with and without timelag correction. From Fig. 6 we see that the 
timelag effect on time domain coherence for lightcurves produced by shots (right panel) is 
much smaller than that for lightcurves sampled from two-dimensional normal distribution 
(left panel). This is because successive bins in a lightcurve consisted of shots have a strong 
connection. 

In summary, we propose a time domain technique in this paper to calculate variability 
coherence for different timescales. In general, the time domain coherence distribution is 
consistent with the Fourier coherence function, but more reliable for the case of low signal 
to noise ratio and more powerful to reveal lower coherence. In practical data analysis, the 
effect of spectral time lags on the time domain coherence can be neglected. This technique 
provides an useful tool for temporal analysis. 

The referee is thanked for helpful comments and suggestions. This research was sup- 
ported by the National Natural Science Foundation of China and made use of data obtained 
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Fig. 1. — Coherence between two simulated lightcurves with a given signal coherence = 
0.64 and different SNR. Prom left to right, SNR = 5, 2 and 1 respectively. Plus sign and 
solid line represent the coherence and its deviation obtained by the time domain technique. 
Circle and dotted line show Fourier coherence function. 
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Fig. 2. — Coherence between two simulated lightcurves with a given SNR = 1 and different 
signal correlation coefficients. Prom left to right, the signal correlation coefficient = 0.81, 
0.64, and 0.36 respectively. Plus sign and solid line represent the coherence and its deviation 
obtained by the time domain technique. Circle and dotted line show Fourier coherence 
function. 
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Fig. 3. — Variability coherence for black hole X-ray binaries. Plus: coherence obtained by 
the time domain technique. Circle: Fourier coherence function. Left panel: Cyg X-1 in the 
low state, PCA/RXTE observation in Dec-17-1996, 2-5 keV vs. 5-13 keV. Middle panel: 
GX 339-4 in the low state, PCA/RXTE observation in Sep-19-1997, 2-6 keV vs. 6-8 keV. 
Right panel: GX 339-4 in the high state, PCA/RXTE observation in Jan-17-1998, 2-5 keV 
vs. 5-7 keV. 
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Fig. 4. — Lightcurve segments of Cyg X-1 in tiie low state, VGA/ RXTE observation in 
Dec-17-1996. Solid line: 2-5 keV, Dashed line: 5-13 keV. Top panel: 10 ms time bin, 
time domain coherence in Fig. 3 is ~ 0.4. Middle panel: 100 ms time bin, coherence ~ 0.9. 
Bottom panel: 1000 ms time bin, coherence is very near to unit. 
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Fig. 5. — X-ray coherence obtained by the time domain technique. Cyg X-1 in the low 
state, PCA/ RXTE observation in Dec-17-1996, 2-5 keV vs. 5-13 keV. Circle: from original 
observed lightcurves. Cross: from lightcurves after spectral timelag correction. Plus: from 
modified lightcurves whose timelag is artificially amplified by 10 times. 
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Fig. 6. — Spectral timelag effect on coherence obtained by the time domain technique. 
Circle : no timelag correction. Plus: after timelag correction. Left panel: coherence between 
two simulated lightcurves based on two-dimensional normal distribution model with timelag 
10ms at all time scale, the given coherence is = 0.64. Right panel: coherence between two 
simulated lightcurves based on shot model with timelag 10ms at all time scale. 



